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Abstract. We study the well-posedness of a model of individual clustering. Given 
p > N > 1 and an initial condition in W 1,P (Q), the local existence and uniqueness of 
a strong solution is proved. We next consider two specific reproduction rates and show 
global existence if N = 1, as well as, the convergence to steady states for one of these 
1 rates. 
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^ ; 1 Introduction 

In [5], a model for the dispersal of individuals with an additional aggregation mecha- 
nism is proposed. More precisely, classical models for the spatial dispersion of biological 
populations read 

d t u = A($(u)) + f(u,t,x). (1) 

where u(t, x) denotes the population density at location x and time t, and f(u, t, x) rep- 
resents the population supply, due to births and deaths. The dispersal of individuals is 
either due to random motion with $(u) = u or rests on the assumption that individuals 
disperse to avoid crowding and $ satisfies 

$(0) = 0, and > 0, for u > 0. (2) 

No aggregation mechanism is present in this model though, as discussed in [5], the onset 
of clustering of individuals in a low density region might balance the death and birth 
rates and guarantee the survival of the colony. To account for such a phenomenon, a 
modification of the population balance (1) is proposed in [5] and reads 

d t u = -V • (u V(u, t, x)) + u E(u, t, x). (3) 

where V is the average velocity of individuals, and E is the net rate of reproduction per 
individual at location x and time t. To complete the model, we must specify how V is 
related to u and E. Following [5], we assume that each individual disperses randomly with 
probability 5 £ (0, 1) and disperses deterministically with an average velocity u so as to 
increase his expected rate of reproduction with probability 1 — 5. The former is accounted 
for by a usual Fickian diffusion ^ while the latter should be in the direction of increasing 
E(u,t,x), say, of the form A X7E(u,t,x) with A > 0. A slightly different choice is made in 
[5] and results in the following system 

d t u = 5 Au — (1 — 5) V • (u u) + u E(u, t, x) , , 

-eAu + uj = XVE(u,t,x). ^ ' 
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After a suitable rescaling, and assuming that the environment is homogeneous, (4) becomes 

dtu = 6 Au — V-(uu) + ru E(u) , 

-eAu + u = E'(u)Vu, ^ 

for x 6 $7 and t > 0, where Q. is an open bounded domain of R N , 1 < N < 3. We 
supplement (5) with no-flux boundary conditions 

n-Vu = n-u = x £ dQ,t>0, (6) 

as suggested in [5]. However, the previous boundary conditions (6) are not sufficient for 
the well-posedness of the elliptic system verified by uj in several space dimensions and we 
must impose the following additional condition given in [3, 4, 8]: 

d n u x n = 0, x £ dQ, t > 0. (7) 

As usual, v x u is the number v\ oj\ + V2 ^2 if N = 2 and the vector field (t^ ^3 — 
V3 u>2, —v\ LO3 + i>3 wi, t>i W2 — V2 (jJi) if = 3. We note that the boundary condition (7) 
is useless if AT" = 1. 



Summarizing, given a sufficiently smooth function E, parameters 5 > 0, e > and 
r > 0, our aim in this paper is to look for (u, u) solving the problem 



5 Au — X7-(uu) + ru E(u), 



(8) 



d t u 

—e Au + u> 

d n u = , u ■ n = 0, 
d n u x n =0, 
u(Q,x) = uo(x), 

In the first part of this paper, we show that, for p > N, the system (8) has a maximal solu- 
tion u in the sense of Definition 2.1 where u e C ([0, T max ), W l *{9.))C\C ((0, T max ), W 2 ' p (n)) , 
see Theorem 2.2. 



x g n,t > 

x € 0,t > 



In the second part, we turn to the global existence issue and focus on space dimension 1, 
and two specific forms of E suggested in [5]: the "bistable case" where E(u) = (l — u)(u—a) 
for some a € (0, 1), see Theorem 2.3, and the "monostable case" E(u) = 1 — u. In both 
cases, we prove the global existence of solution. In addition, in the monostable case, i.e 
E(u) = 1 — u, thanks to the Liapunov functional 

L(u) = J (u log u — u + 1) dx, 

we can study the asymptotic behaviour of solutions for t large, and show that the solution 
u converges, when t goes to 00, to a steady state in L 2 (— 1, 1), see Theorem 2.4. 

In the third part, we investigate the limiting behaviour as e — > 0. Heuristically, when 
e goes to zero, the velocity to becomes sensitive to extremely local fluctuations in E(u), 
and the system (8) reduces to the single equation 

tU = V- {{5-u E\u)) Vu) +r u E(u). (9) 

Clearly (9) is parabolic only if 5 — u E'(u) > for all u > 0. This is in particular the 
case when E(u) = 1 — u, see Theorem 2.6. But this limit is not well-posed in general. 
As a result the population distribution may become discontinuous when neighbouring 
individuals decide to disperse in opposite direction, that is in particular the case when 

E(u) = (1 -«)(«-£»). 
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2 Main results 



Throughout this paper and unless otherwise stated, we assume that 

E £ C 2 (R), 5>0, e>0, r > 0. 
We first define the notion of solution to (8) to be used in this paper. 

Definition 2.1. Let T > 0, p > N , and an initial condition uq £ W 1,P (Q) . A strong 
solution to (8) on [0, T) is a function 

ueC ([0, T), W 1 * (O)) n C ((0, T), W 2 ' P (Q)) , 

such that 

dtu = 5 Au — V ■ (u u u ) +r u E(u), a.e. in [0, T) x £1 
u(Q,x) = uq(x), a.e. in (10) 

d n u = 0, a.e. on [0,T) x 90, 

where, for all t £ [0, T), u u {t) is the unique solution in W 2 ' P (Q) of 

f -eAu u {t) + u u (t) = VE{u(t)) a.e. in fl 

\ u) u {t) ■ n = d n u u (t) x n = a.e. on 30 

Our first result gives the existence and uniqueness of a maximal solution of (8) in the 
sense of Definition 2.1. 

Theorem 2.2. Let p > N and a nonnegative function u$ € W 1,P (Q). Then there is a 
unique maximal solution u G C ([0, T max ), W 1,p (0)) n C ((0, T max ), 1F 2 ' P (0)) to (8) m i/ie 
sense of Definition 2.1 . for some T max £ (0,oo]. In addition, u is nonnegative. 
Moreover, if for each T > 0, £/iere is C(T) smc/i i/iai 

||«(*)Hwi* < for all t £ [0,T] n [0,T max ), 

then T max = oo. 

The proof of the previous theorem relies on a contraction mapping argument. 

We now turn to the global existence issue and focus on the one dimensional case, where 
E(u) has the structure suggested in [5]. 

In the following theorem we give the global existence of solution to (8) in the bistable case, 
that is when E(u) = (1 — u)(u — a), for some a £ (0, 1). 

Theorem 2.3. Assume that u$ is a nonnegative function in W 1,2 (— 1, 1), 

and E(u) = (1 — u)(u — a) for some a £ (0, 1). Then (8) has a global nonnegative solution 

u in the sense of Definition 2.1. 

The proof relies on a suitable cancellation of the coupling terms in the two equations 
which gives an estimate for u in L°°(L 2 ) and for u in L? (W^ 1 ' 2 ). 

Next, we can prove the global existence of a solution to (8) in the monostable case, 
that is, when E{u) = 1 — u, and we show that the solution converges as t — > oo to a steady 
state. More precisely, we have the following theorem 
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Theorem 2.4. Assume that u$ is a nonnegative function in W 1,2 (— 1, 1), 

and E(u) = 1 — u. There exists a global nonnegative solution u of (8) in the sense of 

Definition 2.1 which belongs to L°° ([0, oo); W 1>2 (— 1, 1)) . 

In addition, if r = 0, 

1 f 1 



lim 

t— >oo 



2 



u(t) — — I uq dx 



l 



0. 



and if r > i/ie solution u(t) converges either to or to 1 in L 2 (—l, 1) as i — >• oo. 

In contrast to the bistable case, it does not seem to be possible to begin the global 
existence proof with a L°°(L 2 ) estimate on u. Nevertheless, there is still a cancellation 
between the two equations which actually gives us an L°° (L log L) bound on u and a I? 
bound on d x y/u. 

Remark 2.5. We note that when N = 1, there is a relation between our model when 
E(u) = 1 — u and r = 0, and the following chemorepulsion model studied in [2] 

dtu = 5 d 2 x u + d x (u d x tp), in (0, oo) x (— 1, 1) 

-ed xx ?p + ip = u, in (0,oo) x (-1,1) (12) 

d x u(t,±l) = d x if>(jb,±l)=0 on(0,oo). 

Indeed, define (p = —d x ip, and substitute it into (12) . Then differentiating the second 
equation in (12) we find 

dtu = S d 2 x u — d x (u ip), in (0, oo) x (—1, 1) 

-£dl x ip + (p = -d x u = d x E(u), in (0,oo) x (-1,1) (13) 
d x u(t,±l) = tp(t,±l) = on(0,oo). 

So that u is a solution to our model. 

When E(u) = 1 — u, the limit e — > is formally justified and (8) takes the qualitative 
form of (1) with <£(u) = 5u + ^ u 2 . In this example though, since E' < 0, the individuals 
dispersing so as to maximise E would seek isolation, and there is clearly no mechanism 
capable of producing aggregation of individuals. This observation is actually consistent 
with Remark 2.5. 

Theorem 2.6. Assume that uq is a nonnegative function in W ,2 (— 1, 1), 
and that E(u) = 1 — u. For e > let u e be the global solution to (8) given by 
Theorem 2.4- Then, for all T > 0, 



e^0 

where u is the unique solution to 



T 

lim / \\u £ (t) -u(t)||jj dt = 0, (14) 



o 



d t u = d 2 x (5 u + \ u 2 ) + r u (1-n), x G (-1, 1), t > 0, 

u(0,x) = u (x), z€(-l,l), (15) 

d x u(t,±l) = 0, t>0. 

Since 5 + u > for u > 0, the previous equation (15) is uniformly parabolic and has a 
unique solution u, see [6] for instance. 

The proof of Theorem 2.6 is performed by a compactness method. 
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3 Preliminaries 



We first recall some properties of the following system, 

— e Au> + u = f, in f2, 

u • n = 0, on dQ, (16) 

d n u x n = 0, on dQ, 

where / G (L P (Q)) N and f2 is a bounded open subset of M. N , N = 2, 3. Let us first consider 
weak solutions of (16). For that purpose, we define 

Wx = {v G (H 2 (Q)) N ; v n = 0, and d n v x n = on <9ft} 

and take as the closure of W\ in (i/ 1 (Q)) Jv . 
If / G (L 2 (Q)) N , the weak formulation for (16) is 



£ 



/ Vcj • Vv + / u ■ v = / / • v, for all v G TV 



(17) 



where Vcj • Vv = ^ Vwj • Vfj and u> • v = ^ Wj 

We recall some results about the existence, regularity and uniqueness of solution for 
(17), see [3, 4]. 

Theorem 3.1. // / G (L 2 (Q)) N , (17) has a unique solution in W and there is C = 
C(tt,N) such that 

\H\w <c ||/|| 2 . 



We next consider strong solutions of (16), that is, solutions solving (16) a.e. in Q. In 
this direction the existence and uniqueness of the strong solution to (16) is proved in [8]: 

Theorem 3.2. /// G (L P (Q)) N with 1 < p < 00, (16) has a unique solution in (W 2 ' P (Q)) N 
and 

\\u>\\ W 2, P < ^ ll/llp, (18) 

where K(p) = K(p, O, N). 

In other words, the strong solution has the same regularity as elliptic equations with 
classical boundary conditions. 

We finally recall some functional inequalities: in several places we shall need the fol- 
lowing version of Poincare's inequality 

\\u\\ w i,p < C (||Vu||p + ||u|| ff ), u£W l > p (n) (19) 

with arbitrary p > 1 and q G [l,p]. Also, we will frequently use the Gagliardo-Nirenberg 
inequality 

N_ _ N 

\\u\\ p < C \\u\\ e wh2 M 1 - 9 , with e = _% % , u G W 1 ' 2 ^) (20) 
which holds for all p > 1 satisfying p (N — 2) < 2 N and q G [l,p). 
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4 Local well-posedness 

Throughout this section, we assume that 

E G C 2 QR), and set E(z) = z E{z) for (21) 
Proof of Theorem 2.2. We fix p > N, R > 0, and define for T £ (0, 1) the set 

X R (T):=\ueC([0,T];W 1 ' p (n)), sup |K*)||k^ < R \ , 

{ *6 [0,31 J 

which is a complete metric space for the distance 

d x (u,v) = sup -v(i)|| w i,p ; (u,v) G X R (T) x X R (T). 

te[o,T] 

For u G Xjj(T), and t G [0,T], the embedding of W l > p {VL) in L°°(0) ensures that 
VE(u{t)) G L p (0) so that (16) with / = VE(u) has a unique solution u u G (W 2 '^))^. 
We then define A(u) by 

A«(i, x) = (e* {5 A) u )(x) + / A > [-V • (« w„) + r E(u)} (s, x) ds, (22) 

Jo L - 1 

for (t,x) G [0,T] x $7, where (e* ^ A ^) denotes the semigroup generated in L p (f2) by 5 A 
with homogeneous Neumann boundary conditions. We now aim at showing that A maps 
X R (T) into itself, and is a strict contraction for T small enough. In the following, (Cj)j>i 
and C denote positive constants depending only on Q, 5, r, e, E, p and R. 

• Step 1. A maps X R (T) into itself. 

We first recall that there is C\ > such that 

IN loo < CilMlw 1 ^ (23) 

and 

\\e t(S A) v\\ w i, P < d \\v\\ w i, P , and ||Ve* (<5 A) v\\ p < C x <H H \\v\\ p , (24) 

for all v G VF 1,p (fi). Indeed, (23) follows from the continuous embedding of W l,p (£l) 
in L°°($7) due to p > N while (24) is a consequence of the regularity properties of 
the heat semigroup. 

Consider u G X R (T), and t G [0,T]. It follows from (24) that 



||A«(t)|| 



p 



< d ||«o||p + /"* ||Ve(*- a W A) (u w u )(a) || p 

JO 

+ r / | | e C*— «)C«5 A) s(«)(s)|| p 



i 



< Ci ||n || p + C\ 5 2 / (t - s) z \\u u u (s)\\ p ds 



o 



+ Ci r / ||£(u)(s)|| p ds 



o 



6 



Thanks to (23), we have 

IK*)||oo < Ci \\u(t)\\ m , P <Rd< C 2 . (25) 
Therefore, using elliptic regularity (see Theorem 3.2) and (25), we obtain 

IK(*)II^,p < ^ \\VE(u(t))\\ P < C ||£'||l~(-c 2 ,c 2 ) R<C. (26) 
Using again (25) along with (26) we find 

||Au(t)|| p < CxWuoWp + C [ (t — s) _ ^||«(s)||oo IKMIIp ds 



o 



+ r / \\u(s)\\ p H^Ks))!^ 
JO 

< Ci|KH P + C / (t-a)-^di» + rri2||£;|Ucc ( _ 02|0a) 

JO 

< Ci lluollp + C t* + CT 

< Ci IMI P + C 3 rl (27) 
(recall that T < 1). On another hand, by (24) we have 

||VA«(*)|| P < C x ||V«o|| P + T ||Ve ( *- s) (5 A) V • (u u;„)(s)||p ds 

Jo 

+ r / ||e (i - s)(5 A) V(u #(u))(s)|| p ds 

JO 

< Ci ||Vu ||p + <H / (t-s)~2 \\{Vu-u; u + u V-w„)(a)||p ds 
Jo 

t 



+ r / ||V(S(u))(s)||pds. 
Jo 

Since u G Xr(T), using (25) we can see that 

r ||V(E(u))||p < r \\E'(u) Vu||p < r ||^|U»(-c 3) C a ) ||V«|| P < C 4j 
which gives that 

||VA«(t)|| p < Ci || Vu || p + <H / (*-*)"* I Klloo || Vu||p ds 

Jo 

+ 5"3 / (t - s)-3(||V • w tt ||p |H|oo) ds + Ci T. 
Jo 

Since 

1 1 1 1 oo <Ci \\u u \\ w i, P <C 5 , (28) 
by (26) and (23), we use once more (25) and obtain that 

||VAu(t)|| p < Ci||Vn ||p + C I (t-s)-^ds + C [ (t — s)~^ ds + C^T 

Jo Jo 



< Ci ||Vn || p + C 6 T-2. (29) 



Combining (27) and (29) we get 



sup HAuCt)!!^ < Ci \\u \\ w i,p + C 7 T2 . 
te[o,T] 

Choosing R = 2 C\ | |^o| I an d T G (0, 1) such that 

Ci ||«o||wn, P + C 7 r3 

we obtain that 

sup ||An(t)|| H /i, P < R. 
te[o,T\ 

It follows that A maps Xr(T) into itself. 
• Step 2. We next show that A is a strict contraction for T small enough. 



Let u and v be two functions in Xr(T). Using (24) we have 

||Au(*)-Au(t)|| p < f Ve**"') ( 5 A ) [-u u u + v u v 

Jo 

+ r [ e (t - s)is A) [E(u) - E(v)} 



ds 



o 



ds 



v 



< Ci <5 _ 2 / (t — s)~2 ||— ii u> u + u u>„|| ds 
Jo 



Note that, by (25) and (28), we have 



E(u) - E(v) 



ds. 



| \u U u — V U} v \ \ p < \\u u u — u u v — V u v + u u v \\ p 

< ||u||oo \\uJ u - U) v \\ p + ||u>u||oo Ik-^llp 

< C \\u u - u v \\ p + C \\u - v\\ p , 
and it follows from Theorem 3.2 and (25) that 

\\u u -u v \\ W 2, P < C \\VE(u) - VE(v)\\ p 

< C \\E'{u) Vu - E'(u) Vv - E'(v) Vv + E\u) Vv\\ p 

< C \\E'\\ L oo { _ CM ||Vu-Vi;|| p 
+ C\\E'(v)-E'(u)\\ 00 \\Vv\\ p 

< C \\Vv -Vu\\p + C \\E"\\ Loo{ _ C2>C2) d x (u,v). 

< C 8 d x {u,v). 

Combining (32) and (31) we obtain 



(t-sy 



u uj u + v u) v \\ p ds 



< C T 2 C 8 d x (u, v) + T2 C 9 d x (u, v) 

< ri do d x {u,v). 



Since u and v are bounded by (25), we have 

r \\E{u) - E(v)\\ p < C\\u-v\\ p . 

Then, we get 

It \\E(u) - E(v)\\ p (s) ds < C n Td x (u,v) 
Jo 

Substituting (33) and the above inequality in (30) we conclude that 

||Au(i)-Au(t)|| p <Ci 2 rl d x (u,v). 
Using again (24), we have 

\\VAu(t)-VAv(t)\\ p < f Ve ( '~ s) (<5 A) [-V • (u u u ) +V- (v u v ) 

Jo 



ds 



ds 



Ve (t-s) (8 A) _ ft( v y 

< 5-3 Ci / (t-a)~3 H-V-tucO + V-^cOH 
■/ 

(34) 



V (£(«)-£(«)) (s) ds. 



Since the mapping 



w 1,p (n) x w l ' p (n) — > w 1,p (n) 



is bilinear and continuous due to p > N, we deduce from (26) and (32) that 
S~z C\ || - V • (u (jJ u ) + V • (v (t) v )\\ p < \\u U3 U — v u) v \\ w i, P 

— C \\ u \\w 1 '^ W^u ~ ^vWw 1 '? 

+ C \\oj v \\ w i,p \\u - v\\ w i, P 

< C d x (u,v) + C \\u — v\\ w x, P 

< C13 d x (u,v) 

Thus, 

Ci / (t -s)~3 II -V- (« + V- (« w„)||p ds < C 13 d x (u,v). 
Jo 

On the other hand, due to (25) and the embedding of W 1,P (Q) in L°°(Q) 

\\V(E(u)-E(v))\\ p = \\E'(u)Vu-E'{v)Vv\\ p 

< C ||i'(^||oo || Vu - V«|| p + C \\E'{u) - £»||oo \\Vv\\ p , 

< C \\E'\\ Laoi _c 2 ,c 2 ) \\Vu-Vv\\ p + C \\E"\\ La o { _ C2 ,c 2 ) lk-Hloc 

< C \\u — v\\ w i, P . 



Then 



r Ci ||V(£(u) - E(v))\\ p < C14 ||n - v\\ w i, P . 
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Therefore, 

\\VAu(t) -VAv(t)\\ p <T$ Ci 3 d x (u,v) +T C u d x (u,v) < C 15 d x (u,v). 
Finally we get 

d x {Au(t) -Av(t)) < Ci 6 d x (u,v). 

Choosing T G (0, 1) such that Ta Ciq < 1 we obtain that A is indeed a strict con- 
traction in Xr(T) and thus has a unique fixed point u. 

Furthermore, since u G C ([0, T], W^ p {n)) andp > N, we have VE(u) G C ([0, T], L p ($7)) 
and we infer from Theorem 3,2 that u M 6 C ([0, T],W 2 > P (Q)) . Comb ining this prop- 
erty with the fact that W l ' p {£l) is an algebra, we realize that both V • (u u u ) and 
u E(u) belong to C ([0, T], L p (£l)). Classical regularity properties of the heat equa- 
tion then guarantee that ueC ((0, T], W 2 ' P (Q)) and is a strong solution to (10). 

• Step 3. Thanks to the analysis performed in Steps 1 and 2, the existence and 
uniqueness of a maximal solution follows by classical argument, see [1] for instance. 

• Step 4. Since clearly solves (10), and no > 0, the positivity of u follows from the 
comparison principle. 

□ 



5 Global existence 

From now on we choose N = 1, O = (—1, 1), p = 2 and we set ip = u: u to simplify the 
notation. 

5.1 The bistable case: E{u) = (1 — u)(u — a). 
In this case, the system (8) now reads 



(35) 



dtu = 5 d 2 x u — d x (u tp) + r u (u — a)(l — u), x G (— l,l),t>0 

-edl x Lp + Lp = {-2u + (a + 1)) d x u, x G (-1, l),t > 

d x u(t,±l) = (p(t,±l) = 0, t>0, 

k u(0,x) = u (x), x G (-1, 1), 

for a some a G (0, 1). 

bince E G C 2 (R), Theorem 2.2 ensures that there is a maximal solution 
u of (35) in C ([0, T max ), W^^-l, 1))) n C ((0, T max ), W 2 >f (-1, 1)) . 

To prove Theorem 2.3 we show that, for all T > and t G [0, T] n [0,T max ), u(t) is 
bounded in W 1,2 (-l,l). 

We begin the proof by the following lemmas which give some estimates on u and <p. 

Lemma 5.1. Let the same assumptions as that of Theorem 2.3 hold, and u be the non- 
negative maximal solution of (35). Then for allT > there exists C\{T), such that u and 
ip satisfy the following estimates 

\\u(t)\\ 2 <d(T), for all t G [0, T] Pi [0, T max ), (36) 
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f \\d x u\\ 2 2 dt < Cx(T) for all iG [0, T] n [0, T max ), (37) 

JO 

and 

/ ||^(t)||^i, a dt < Ci(T) for all i € [0,T] n [0,T max ). (38) 
■/ o 

Proof. Multiplying the first equation in (35) by u(t) and integrating it over (—1, 1), we 
obtain 



d_ 

dt 



J \u\ 2 dx = -2sj \d x u\ 2 dx + 2 J u(pd x udx + 2r J u 2 E{u) dx. (39) 

Multiplying now the second equation in (35) by ip and integrating it over (—1,1) we obtain 

e J \d x <p\ 2 dx + j \<p\ 2 dx = —2 J u ip d x u dx + (a + 1) J d x u <p dx. (40) 

At this point we notice that the cubic terms on the right hand side of (39) and (40) cancel 
one with the other, and summing (40) and (39) we obtain 



d. 
d7 



\u\\l+e\\d x cp\\l + \\cp\\l + 2 5\\d x u\\l = 2r J u 2 E(u) dx + (a + l) J d x 



u (p dx. (41) 



We integrate by parts and use Cauchy-Schwarz inequality to obtain 

(a + 1) / d x u ip dx = — (a + 1) / u d x ip dx < - — ^ — — ||n||2 + - H^WHl- 
J-i J -i 2 £ 2 

On the other hand, it 2 E(u) < if u ^ (a, 1) so that 



J u 2 E(u) dx < 2 (1 - a) 



The previous inequalities give that 

4|Ml! + \ ll^lll + IMIl + 2 S \\d x u\\ 2 < \\u\\l + 4 r (1 - a), 

at 2 2 e 

Therefore, for all T > there exists C X (T) such that (36), (37) and (38) hold. □ 

Lemma 5.2. Let the same assumptions as that of Theorem 2.3 hold, and u be the non- 
negative maximal strong solution of (35). For all T > 0, there is Coo(T) such that 

IK*)lloo < Coo(T), for all t G [0,T] D [0,T max ). (42) 

Proof. The estimates (36) and (37) and the Gagliardo-Nirenberg inequality (20) yield that 
there exists C2 (T) such that 

ft ft x 3 

/ \\u d x u\\ 2 ds < C / ||n|| 2 2 H^ulll (is < C 2 (T) for all t G [0,T] n [0,T max ). 

The second equation in (35) and classical elliptic regularity theory ensure that there exists 
C (T) such that 

f \\<p\\w2,3 ds < C{T) for all t G [0, T] n [0,T max ), 

JO 
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which gives in particular, since W 2,2 (—l, 1) is embedded in W 1), 

f | \dM |oo ds < C 3 (T) for all t G [0, T] n [0, T max ). (43) 

JO 

Now, we multiply the first equation in (35) by q u^ 1 where q > 1 and integrate it over 
(—1, 1) to obtain 



^ C \u\ q dx = -4^ — ^ f 1 \d xU 2\ 2 dx + (q-l) [\d x u q 



+ r g / u 9 (u — a) (1 — u) dx 



i 



< -(«-!)£ 



dx + r g 2 (1 — a) 



Using Holder's inequality , we obtain 

_d 
di 

Introducing 



_|| u ||9<(g-l) p^lloo ||u||« + 2gr. (44) 



<f>(t) = [ P^(s)||oo ds < C 3 (T) for all t G [0,T] D [0,T max ), 
JO 

the bound being a sequence of (43), we integrate (44) and find 



\\u(t)\\ q g < \\u \\ q e^" 1 ) +2 q r / 

■/ o 

< (\\uo\\ q + 2qr)Te qC ^ T \ 
\\u(t)\\ q < ((\\u \\ q + 2qr) T)< e c ( T ) for all i G [0, T] n [0, T max ). 
Consequently, by letting g tend to oo, we see that there exists Coo(T) such that 
IK*) I loo < Coo(r), for all t G [0,T] n [0,T max ). 



□ 



Lemma 5.3. Let £/ie same assumptions as that of Theorem 2.3 hold, and u be the non- 
negative maximal strong solution of (35). For all T > 0, there is C^(T) such that 

\\d x u(t)\\ 2 < C A (T) for all t G [0,T] n [0,T max ). (45) 

Proof. We multiply the first equation in (35) by (—d xx u) and integrate it over (—1, 1) to 
obtain 



1 

2 dt 



J \d x u\ 2 dx = —5 J \d xx u\ 2 dx + J d x (u tp) d xx u dx 

+ r j u (1 — u) (u — a) {—d 2 x u) dx 

= —5 J \d xx u\ 2 dx + J (u d x ip + d x u tp) d xx u dx 

+ r J u (—u 3 + (a + l)u 2 — a) d xx u dx. 
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1 d 

2 di 



Using Cauchy-Schwarz inequality and Lemma 5.2 we obtain, 

J \d x u\ 2 dx < —5 J \d 2 x u\ 2 dx + - J \d 2 x u\ 2 dx + C WuWf^ J \d x ip\ 2 dx 

+ ljjd 2 xx u\ 2 dx + C\Ml J\d x u\ 2 dx 
5 I' 1 

+ n \dl x u\ 2 dx + C \\{u (-v? + (a + l)u 2 - a) 
3 7-i 



oo 



< C(T) J \d x tp\ 2 dx + C \\ip{t)\\lo f \d x u\ 2 dx + C(T). (46) 



□ 



Using (38) and Sobolev embedding theorem we obtain the following estimate 

/" IMIL ds < C(T) for all t G [0,T] n [0,T max ). (47) 
7 o 

Since (38) and (47) hold, then it follows from (46) after integration that 
\\d x u(t)\\ 2 < C 4 (T), for all t € [0, T] n [0,T max ). 

It remains to prove Theorem 2.3. 

Proof of Theorem 2.3. For all T > 0, Lemma 5.3 and the estimate (36) ensure that 

\\u{t)\\ w i,2 < C(T), for all t G [0,T] n [0,T max ), 

which guarantees that u cannot explode in W l,2 {— 1, 1) in finite time and thus 

that T max = oo. □ 

5.2 The monostable case: E(u) = 1 — u 
For this choice of E, system (8) now reads 

dtu = 5 d xx u — d x (u (p) + r u (1 — u) x G (—1, 1), t > 

-e^y + V = - §b«i xG(-l,l), t>0 , . 

d x u(t,±l) = (p(t,±l) =0 t>0, 1 J 

, u(0,x) = n (x) x G (-1, 1), 

bmce E G C 2 (R), Theorem 2.2 ensures that there is a maximal solution u of (48) in 
C ([0, T max ),W^(-l, 1)) n C ((0,T max ), TU 2 'f(-l, 1)) . 

In contrast to the previous case, it does not seem to be possible to begin the global 
existence proof with an L°°(L 2 ) estimate on u. Nevertheless, there is still a cancellation 
between the two equations which actually gives us an L°°(Llog-L) bound on u and a L 2 
bound on d x y/u. Integrating (48) over (0,t) x (—1,1) and using the nonnegativity of u, 
we first observe that, 

IKi)||l < |K||i + 2r t, for all i G [0,T max ). (49) 
To prove Theorem 2.4 we need to prove the following lemmas: 
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Lemma 5.4. Let the same assumptions as that of Theorem 2.4 hold, and let u be the 
maximal solution of (48). Then for all T > 0, there exists a constant C\(T) such that 

[ \\d x ^\\ 2 2 ds < Ci(T), for all t G [0,T] n [0,T max ), (50) 
J o 

and 

f \\<p\\ 2 w i,2 ds < d(T), for all t € [0,7] n [0,T max ). (51) 
J o 

Proof. The proof goes as follows. On the one hand, we multiply the first equation in (48) 

by (log it + 1) and integrate it over (—1, 1). Since u (1 — u) logu < and u (1 — u) < 1, 

d f 1 f 1 1 f 1 

— / u logn dx = — (5 d x u — u tp) (— c^tt) fix + r / u (1 — u) (logn + 1) c£c 
di y_i J -i u J_ x 

< - f - (d x u) 2 dx+ [ if d x u dx + 2 r. (52) 

7-i « J -i 

On the other hand, we multiply the second equation in (48) by ip and integrate it over 
(— 1, 1) to obtain 

e J \d x (p\ 2 dx + J \<p\ 2 dx = — J d x u ip dx. (53) 
Adding (53) and (52) yields 

u logu dx + £ 115^111 + ll^li < -4 8 J \d x y/u\ 2 dx + 2 r. (54) 

Finally, (50) and (51) are obtained by a time integration of (54). □ 

Lemma 5.5. Let the same assumptions as that of Theorem 2.4 hold, and let u be the 
maximal solution of (48). Then for all T > 0, there exists a constant (^(T) such that 

\\u(t)\\ 2 <C 2 (T), for all t £ [0,T] n [0,T max ), (55) 

and ^ 

[ \\d x u\\ 2 2 ds < C 2 (T), for all t € [0,T] n [0,T max ). (56) 
Jo 

Proof. A simple computation shows that, since u 2 (1 — u) < 1, 

J u 2 dx < —2 5 j \d x u\ 2 dx + 2 j u ip d x u dx + 4 r. (57) 



d 
ch 7 



Using Cauchy-Schwarz inequality, Gagliardo-Nirenberg inequality (20), Young inequality 
and (49) we obtain that for all T > 0, 

2 j up d x u dx = — j u 2 d x p dx < \\u\\ 2 \\d x <fi\\2 

< c ( |H|^ li2 \\ u \\f) \\d x p\\ 2 < c(T) \\d x p\\ 2 \\u\\ w i,2 

iw 1 



< 5 |H|Li, 2 + C(T) \\d x p\\ 2 
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We substitute the previous inequality in (57) to obtain 

^\\u\\ 2 2 <-2 5\\d x u\\ 2 2 + 5\\u\\ 2 wl , 2 + C(T) \\d x ip\\l + Ar. (58) 

Integrating (58) in time, and using (51) yield that there exists C^{T) such that (55) and 
(56) hold. □ 

Now we are in a position to show the global existence of solution to (48). 

Proof of Theorem 2.4 (global existence). 

By elliptic regularity, and the continuous embedding of W 2,2 (— 1,1) in W 1,co (— 1, 1), 
we have 

P^(t)||oo < C Mt)\\ W 2a < C \\E'{u) d x u\\ 2 < C \\d x u\\ 2 , 
which together with (56), implies that 

f\\dM.s)\\lo ds<C f \\d x u{s)\\l ds<C(T), for all t G [0,T] n [0,T max ). 
Jo Jo 

Thanks to this estimate, we now argue as in the proof of Lemma 5.2 and Lemma 5.3 to 
get that 

IHOHoc + \\d x u(t)\\ 2 < C(T), for all t G [0,T] n [0,T max ). 
Thus, the maximal solution u of (48) cannot explode in finite time. □ 

To complete the proof of Theorem 2.4, it remains to prove the asymptotic behaviour 
of u when t — > oo. We note that we have the following lemma which controls the L l {— 1, 1) 
norm of u. For / G L (— 1, 1), we set 



1 



2 



< / >= n / /0) dx 



-1 



Lemma 5.6. Let the same assumptions as that of Theorem 2.4 hold, and let u be the 
nonnegative global solution of (48) . For r > 0, there exists a constant Co > such that 

<< u(t) >< C , t G (0, oo), (59) 

and if r = 

< u(t) > = < u >= -\\u \\i, tG(0,oo). (60) 
Proof. We note that if r = 0, 4: < u >= 0, so that 

< u (t) >= 2 IK*)lli = \ INIIi- 

If r > 0, 

— < u(t) > = r < u(t) > -r < u 2 (t) > 
< r < u(t) > -r < u(t) > 2 , 

whence < u(t) >< max{l,< uq >}. □ 

Next we turn to the existence of a Liapunov functional for (48) which is the cornerstone 
of our analysis. 
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Lemma 5.7. Let the same assumptions as of that Theorem 2.4 hold, and let u be the 
nonnegative global solution of (48). There exists a constant C\ such that 



OO fl 



u (u-1) logudx dt+ / ( 113^11! + e \\d x <p\\l + M\l ) dt < C u (61) 
/[) J-l Jo 

and 

poo 

~,ds<C x . (62) 



Proof. Let us define the following functional L 



L(u) = j (u log u — u + 1) dx > 0, 



and show that it is a Liapunov functional. Indeed 
d 



f 1 \dxu\ 2 f 1 f 1 

/ — dx + / (p d x u dx + r I u logu (1 — u) dx. (63) 

J-i u J_ t 7_i 



Combining (63) and (53) we obtain that 

^ L(u) = —4 5 J \d x y/u\ 2 dx — e J \d x ip\ 2 dx — J \ip\ 2 dx + r J u logu (1 — u) dx. 



dt 



= -D(u, if) < 0, (64) 
since u logu (1 — u) < 0. Then for all t > 

£(«(*)) + / D(«(a),y>(a)) < L(n ). 
■/ o 

Since an d ^ are nonnegative , we have 



L(u ) < J (u 2 + 1) dx < \\u \\l + 2, 



and 

2 f 1 

L(u(t)) > / uq dx, 

e J-i 

so that 

/"* 2 

/ D(u(s), <p(s)) cis < 1 + lluollo H h2 < u >, t>0. 

Jo e 

Therefore, (65) yields there exists C\ > such that 

/"OO 

/ L>(u,9?) (it < Ci. 



(65) 



(66) 



From (66), we see that (61) holds true. In addition, inequality (61) together with Sobolev's 
embedding theorem give (62). □ 

In the following lemma we show that {u(t) : t > 0} is bounded in W ' 2 (— 1, 1). 

Lemma 5.8. Let the same assumptions as that of Theorem 2.4 hold, and let u be the 
nonnegative global solution of (48). Then u belongs to L°° ((0, oo); W 1,2 (— 1, 1)). 
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Proof. First, 



= 1 1^1 ll — 4 r < u > 2 +4 r < u > < u 2 > 
di 

< ^ IMll + 2r C |H||. (67) 

Multiplying the first equation in (48) by 2 u, integrating it over (—1,1), and using the 
Cauchy-Schwarz inequality and the fact that u > we obtain 

"'"III — —2 5 I \d x u\ 2 dx — I u 2 d x ip dx + 2 r I u 2 (1 — u) dx 



. /IIM . < -2 5 J \d x u\ 2 dx-J u 2 d x tpdx + 2rj 

< -2 8\\d x u\\l + \\u\\l\\d x ip\\ 2 + 2r. (68) 



Gagliardo-Nirenberg inequality (20) together with the Poincare inequality (19) and (59) 
give 

IMI4 < C \\u\\ w i,2 \\u\\i < C ( \\d x u\\ 2 + \\u\\i ) 

< C(\\d x u\\ 2 + 1). (69) 

Thus 

IMI2 \\dM\2 < c (\\d xU \\ 2 + 1) ||a^|| 2 . (70) 

Substituting (70), (69) and (68) into (67), and using Young and Holder inequalities to 
obtain 

^||«-<«>||1 < -26 \\d x u\\ 2 2 + C {\\d x u\\ 2 + \) ||a^|| 2 + 2r 

2 4 
+ 2r C \\u\\l 1 1 

< -2 6 \\d x u\\ 2 2 + C \\d x y\\ 2 2 + 6 \\d x u\\ 2 2 + C. 
Using Poincare's inequality we get 

^ \\u(t)- < u{t) > + a \\u{t)- <u(t) > Hi < C\\dMt)\\l + C, 
for some a > independent of t. Integrating this differential inequality gives 

<«(*)> Hi < \\uo-<u >\\ 2 2 + [ e as (C\\dMs)\\ 2 + C) ds 

Jo 

\\u(t)- <u(t) > 111 < Ce~ at + C [ e a (s -V (\\d x <p{s)\\ 2 2 + 1) ds 

Jo 

Since e~ a * < 1, and e a < 1 as s < t we obtain 

rt 1 
\\u(t)- <u(t) > Hi < C + C \\d x (p(s)\\ 2 2 ds + -, 

Jo a 

Using (61) we end up with 

\\u(t)- < u(t) > \\l < C, (71) 
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where C is independent of t. Therefore, u belongs to L°°((0, oo); L 2 (— 1, 1)). 

It remains to show that d x u is in L°°((0, oo); L 2 (— 1, 1)). We multiply the first equation 
in (48) by —d xx u and integrate it over (—1,1). Since u > we use Cauchy-Schwarz and 
Young inequalities and (71) to obtain 



Id 

2dt 



H^a^Hl < — 8 ||o*zx u lli + J dl x u (d x u (p + d x if u) — r J u (1 — u) d xx u d.r 
+ cj \u\ 2 \d x ip\ 2 dx - r J ud xx udx-2r J u\d x i 



r u\ 2 dx 



< -5 \\d 2 xx u\\ 2 + \ \\d x y\\ 2 \\d x u\\l + S - \\d 2 xx u\\ 2 

+ CIMIL \\d x v\\l + C+~ \\d 2 xx u\$. (72) 

1 2 

Since d x u G W^' (-1,1), using the Gagliardo-Nirenberg inequality (20) and the classical 
Poincare inequality (19). we obtain 

\\d x u\U<C\\d 2 xx u\\l \\d x u\\l (73) 



Then, we substitute (73) into (72), and by Young inequality, the Sobolev embedding, (59) 
and (71) . we obtain 

\^ t \\d x u\\ 2 2 < ~ ll^«Hl + C||Ml2 IIC«I|2 ||9««||x + C|H|^ ||^||i + C 

< -\ \\d 2 xx u\\l + c\\dM\ 2 \\d xU \\ 2 + c\\ u \\ 2 wl , 2 \\d xV \\ 2 + c 

< -~ \\d 2 xx u\\l + C \\d x u\\l \\dM\ 2 + C (P^lli + 1). (74) 
Since d x u G W M (-1,1), we use once more the classical Poincare inequality to obtain 

±\\d x u\\ 2 +f3\\d x u\\ 2 < C \\d x u\\ 2 \\dM\ 2 + C [\\dM\ 2 + 1), 
for some (3 > independent of t. 
Define 

cf>(t) = [ \\dMs)\\ 2 ds, t>0, 
Jo 

and notice that , since Hc^H 2 . belongs to L 1 (0,oo) by (61), 

poo 

0<<t>(t) <<Poo= / \\dMs)\\ 2 ds. 
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Integrating the previous differential inequality we find 

\\d x u(t)\\ 2 2 < ll^uolll e^W-** 

+ C7 f (1 + d>'(s)) ef*-t)+c<Kt)-CK.) ds 



< \\d x u \\ 2 2 e c ^+C e c ^ f\e^ s ^ + 0'(s)] 

< H^uolll e c *nCe c ^ fin- 



ds 



Therefore, c^u belongs to L°° ((0, oo), L 2 (— 1, 1)), and Lemma 5.8 is proved. □ 

Lemma 5.9. Let the same assumptions as that of Theorem 2.4 hold, and let u be the 
nonnegative global solution of (48). There is C*2 such that 

poo 

/ \\d t u\\ 2 dt<C 2 . (75) 
J o 

Proof. We multiply the first equation in (48) by dtu and integrate it over (—1, 1) to obtain 
J (dtu) 2 dx = 5 J d 2 x u dtu dx — J d x (u (p) dfU dx + r J u (1 — u) dtu dx 



2 ^ dxU ^~ J (dxU (p + u d x ip) d t u dx + r J 



(u — u 2 ) dtu dx. 



Using Young and Cauchy-Schwarz inequalities we obtain 



\\d t u\\l < J-±.\\d x u\\l + J^\ dx u\ 2 \ip\ 2 dx + ^\\d t u\\l + J^\u\ 2 \dM 2 dx + )\\d t u\\ 2 



dt 

which gives 



d f 1 (u 2 u 3 \ , 
+ r — — — ) dx, 



di [ 2 



\d x u\\ 2 + J F(u) dx} + \\\d t u\\ 2 < \\d x u\\ 2 |M|L + \\u\H \\d x <p\\l 



where F(u) = r {-^- + ^) > -§. 

Next we integrate the above inequality in time, and use (62), (61) and Lemma 5.8 to 
obtain 



\\d t u\\ 2 ds<C + C [ (IMlL + IIMl!) ds<C 
Jo 

for t > where C is independent of i. We have thus proved (75). □ 



r 1 
3 + 2 



To end the proof of Theorem 2.4, our aim now is to look at the large time behaviour 
of the solution. 

Proof of Theorem 2.4, (large time behaviour). In this proof, we follow [7]. 
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By Lemma 5.8, the family {u(t), t > 0} is bounded in W ,2 (— 1, 1)- Since the embed- 
ding of W 1 > 2 (— 1,1) in L 2 (— 1,1) is compact then, there are a sequence of positive time 
(t n ), such that t n — > oo, and z S L 2 (— 1, 1) such that 



Consider 
and 

We first prove that 



z = lim u(t n ) in L 2 ( — 1, 1) and a.e. in (—1, 1). 
C/ n (s, x) = u(i n + s, x), x G (—1, 1), — 1 < s < 1, n > 0, 

$n(s,x) = <p(t n + S,x), -1 < S < 1. 



U n — >z as n -)• oo, in C ([-1, 1]; L 2 (-l, 1)) . 
Indeed for each s G (—1, 1) 



(76) 



1 /-tn + l 



|it(t n + s, x) — u(t n , x)| 2 dx < 



\d t u\ 2 dt dx. 



i Jt n -i 



Hence 



sup \\U n (s) - u(t n )\\ 2 < 

se[-i,i] 



1 fOO 




\dtu\ 2 dt dx 



-l Jt n -i 

The right hand side goes to zero as n — > oo by Lemma 5.9. Letting n — > oo in the above 
inequality gives (76). 

Next, using the definition of D(u, ip) which is given in (64) we obtain that 
J ^(r\\U n log U n (U n - 1)111 + 11^^111 + ||$n||| + e H^nlll) ds 



< 



tn+l 

t„-l 
oo 



< 2 



r \\u(s) logu(s) ( u (s)- 1)11, + \\d xV ^{s)\\ 2 + M S )\\ 2 + e \\d x <p(t 



D(u, (p) ds. 



ds 
(77) 



t„-i 



The right-hand side of (77) goes to zero as n — > oo by (66), so that 

$ n — >0 asra^oo, in L 2 ((—1, 1); W 1,2 (—l, 1)) . 
In addition, using Cauchy-Schwarz inequality, (59) and (77) we obtain 



\\d x U n (s)\\l ds 



1 / rl 



y/U n {s) \d x y/U n {s)\ dx) ds 



< 4 



l V-i 
l 



ll^nWHl \\dxVUn{s)\\l ds 

1 

/•oo 

< C D(u,ip) ds. 
Jt n -i 

Since the right-hand side goes to zero as n — > oo by (66), then we have 
d x U n — >0 asrwoo, in L 2 ((— 1, 1); L 1 (— 1, 1)) . 



(78) 
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Since the limit in the sense of distribution is unique, (76) and (78) yield that 

8 x z = 0. (79) 

If r = 0, (79) together with (60) and (76) give that z =< uq >. We have thus shown 
that < no > is the only cluster point of {u(t), t > 0}. Since {u(t), t > 0} is relatively 
compact in L 2 (— 1, 1) thanks to its boundedness in W 1,2 {— 1, 1) (see Lemma 5.8), we con- 
clude that u(t) converges to < uq > in L 2 (— 1, 1) as t — > oo. 



If r > 0, by (77) we have 

J \\U n logU n (I7 n -l)||i ds^O, asn^oo, (80) 

Since (U n ) is bounded in L°° ((—1, 1) x (—1, 1)) thanks to the boundness of {u(t), t > 0} 
in W^ 2 (-l, 1) and the embedding of W l ' 2 (-l, 1) in L°°(-l, 1), we infer from (76), (79), 
(80) that z \ogz (z — 1) = 0, that is z = or z = 1. Therefore and 1 are the only two 
cluster points of {u(t), t > 0} as t — > oo. Since the w-limit set of u is a compact connected 
subset of L 2 {— 1,1), see [1, Theorem 9.1.8] for instance, we conclude that u(t) converges 
either to or to 1 in L 2 (—l, 1) as t — > oo. 

□ 

6 Limiting behaviour as e — > 

When E(u) = 1 — u, letting e — > in (48) formally leads to (9) which is well-posed 
since E' < and 5 > 0. The purpose of this section is to justify rigorously this fact and 
prove Theorem 2.6 . Let T>0,5>0,r>0,e>0 and a nonnegative initial condition 
no G W ,2 (— 1, 1). We discuss the limit as e — > of the unique solution u £ of 

d t u £ = 6 d 2 x u £ - d x (u £ ip £ ) + r u £ (1 - u £ ) in (0, T) x (-1, 1), 

u £ (0,x) = u (x) in (-1,1), (81) 

d x u £ (t,±l) = on(0,T), 

given by Theorem 2.4, where (p £ is the unique solution of 

-e d 2 xx ip £ + ip £ = -d x u £ in (0, T) x (-1, 1), 
<p E (t,±l) = on(0,T). 

6.1 Estimates 

Lemma 6.1. There is C\{T) independent of e such that 

[ {5 WdvyfeWl + £ 11^^111 + ||^ £ |||) dt< Cx(T), (82) 
Jo 

Proof. By (54) see (the proof of Lemma 5.4), we have 



d_ 

dt 



J u £ \ogu £ dx + e \\d x ip £ \\\ + \\tp £ \\l < -4 5 J \d x ^/u~ £ \ 2 dx + 2 r, 



from which (82) follows by a time integration. □ 
Using Gagliardo-Nirenberg inequality (20) we obtain the following estimate: 
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Lemma 6.2. For 2 < p < 6, there exists C2(T,p) independent of e such that 

| Kill dt < C 2 (T,p). (83) 



Proof. For t G (0,T), thanks to (82), we can use Gagliardo-Nirenberg inequality (20) on 
y/v^ and we obtain for all p G [2, oo) 



where 



Therefore 



Hv^IIp < c \\V^S)\\ e w ^ \\V^it)\\\- e , 

2 p 



(84) 



p-2 



WV^WM < c||V^)lUi, 2 IIV^OOH 



p+2 



12 
p+2 



Since 



K(t)||I < ciiv^llJ^ IkWHi 4 • 



M*)lli < IKHi + 2 rT, 



and < 2 for 2 < p < 6 , the estimate (83) follows from (82) and the previous 
inequalities. □ 

Lemma 6.3. There is C 3 (T) independent of e such that 



T 

\\d x u e \\\ dt<C 3 (T). 



Proof. Holder and Young inequalities together with (82) and (83) with p = 6 yield 



(85) 



\\d x u e \\l dt < 2 

2 



u £ d x y/u^\\l dt 

2 

T 3 3 

ll^eHg 1 1 d x yju e 1 1 2 dt 



< 22 

< C (IKHI + 115^111) dt<C 3 (T), 

Jo 

which gives the result. 

Lemma 6.4. There is C±(T) independent of e such that 



□ 



J^\\d t u £ \\ {w2 3 y dt<C A (T). 



Proof. Consider ip £ W 2 '2(-l, 1) and t £ (0,T), we have 



dtu e V dx 



J [d x (5 d x u e - u e ip e ) + r u e E{u £ )} ip dx 



(-5 d x u e d x ip + u e <p s d x ip + r u e (1 - u e ) if>) dx 



-i 



< S H^VHoo H^Uelll + Halloo 1 1 "tie- 1 1 2 \\fe\\2 + r \\u e (l - U, 



e)\\l II ^ II 00 
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Using the embedding of W '2 (—1, 1) in W 1,0 °(—l, 1), and Young's inequality, we end up 
with 



J d t u £ if) dx < (5 \\d x u £ \\i + \\u £ \\ 2 \\feW2 + r \\u £ \\i + r \\u £ \\f) 
< C (\\d x u £ \\3 + \\u e \\% + \WeW2 + 1. 
and a duality argument gives 



ll^(*)ll (wa ,t y < C (J|cU £ ||| + ||u e ||£ + \\<p e \\i + 1 
Integrating the above inequality over (0, T) and using Young's inequality we obtain 

J \\d t u £ (t)\\ {w ^s y dt<C(T)^ (\\d x u e \\l + \\u £ \\l + \\<p £ \\l + lj dt. 

By Lemma 6.1, Lemma 6.3 and Lemma 6.2 with p = 4 the right-hand side of the above 
inequality is bounded independently of e and the proof of Lemma 6.4 is complete. □ 

6.2 Convergence 

In this section we discuss the limit of u £ as e — > 0. For that purpose, we study the 
compactness properties of (u e ,<p £ ). 

Proof of Theorem 2.6. Thanks to Lemma 6.2 and Lemma 6.3, (u £ ) £ is bounded in 
Ll((0,T);^ 1 't(-l,l)) while (d t u £ ) £ is bounded in L 1 ((0,T); (W 2 'l)'(— 1, 1)) 
by Lemma 6.4. Since W 1,1) is compactly embedded in C([— 1,1]) and C([— 1,1]) 



is continuously embedded in (W 2, 2 )'(— 1, 1), it follows from [9, Corollary 4] that (u, 
is relatively compact in L2 ((0, T); C{[— 1, 1])) . Therefore, there are a sequence (ej) of 
positive real numbers, £j — > 0, and u £ ( (0,T); W 1 '^ ) such that 



eye 



3 / , , 1 3 



u £j inL2 ^(0,r);W i '5(-l > l)J , (86) 

and 

u £j . ->• u in L§ ((0,T);C[-1,1]) and a.e. in (0, T) x (-1, 1). (87) 
Since (u £ ) £ is bounded in L°° ((0, T); L 1 (— 1, 1)) by (59), it follows from (87) 

rT r T 3 a rT 



3 
2 

00 



/ /" M- 11- r 

/ - n|| 2 (it < / ||u £j . - u\\l \\u £j - u\\Zo dt < 

Jo Jo Jo 

when £j — > 0. In particular, we have 

u £j — >«, inL 2 ((0,T) x (-1,1)). 
Observe that the nonnegativity of u follows easily from that of {u £ . ) by 

Owing to Lemma 6.1. we may also assume that 

ip £j ip in L 2 ((0,T) x (-1, 1)) as Ej 0, (89) 
£j d x ip Sj -> in L 2 ((0, T) x (-1, 1)) as Ej -»■ 0. (90) 
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It remains to identify the equations solved by the limit u of (u £ ). Let ip G C 2 ([0,T] x 
[-1,1]) with ip(T) = 0. Since 

T [-1 t-T [-1 

/ d t u £j ip dxdt = / ((-5 d x u £j + u £j cp ej ) d x ip + r u £j E(u £j ) ip) dxdt (91) 
o J-i Jo J-i 

and 

rT rl rT rl rT rl 

£j / / d x ip £j d x if> dxdt + / ip £j ip dxdt = - I / d x u £j -0 dxdt. (92) 
Jo J-i Jo J-i Jo J-i 

Owing to (86), (89) and (90), it is straightforward to pass to the limit as Ej — > in (92) 
and find 

J ip ip dxdt = — / d x u ip dxdt, 
'0 J-i Jo J-i 

which gives that 

ip = -d x u. (93) 

Next, by (87) and (86) we see that 

r-T r l rT rl 



J-l 



dtu £j ^ dxdt — > — I u df>P dxdt — \ uq(x) tp(0,x) dx, 
Jo J-i J-i 



as Ej — > 0, 



and 

r-T rl rT r l 



/ / d x u £j d x %p dxdt — > j j d x u d x ip dxdt as Ej — > 0. 
Jo J-i Jo J-i 



From (88), (89) and (93) we see that 



T rl rT rl 

u £i <Pei 9 x tp dxdt — 



J rc 3 

J-l 



u d x u d x tp dxdt, as Ej — > 0. 

Jo J-i 



From (88) we get 

cT rl rT r l 



r 



I I u £j E(u £j ) ip dxdt — > r I I u E(u) ip dxdt, as Ej — > 0. 
Jo J-i Jo J-l 



Thus we conclude that u satisfies 



T ,T rl 

(dtu, ip) dt = / ((—5 d x u — u d x u) d x ip + r u E(u) VO dx dt, 

for all test functions ip. Therefore, u is a weak solution of (15), and classical regularity 
results ensure that u is actually a classical solution of (15). Since it is unique and the only 
possible cluster point of (u £ ) £ in I? ((0,T) x (—1, 1)), we conclude that the whole family 
(u £ ) £ converges to u in L 2 ((0,T) x (—1, 1)) as e — > 0. □ 
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